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Abstract.  Different  methods  of  calculation  of  clusters  are  concerned  from  thequantum  chemistry  point  of  view.  Results 
obtained  with  the  help  of  these  methods  for  carbon,  water  sulfur, alkane  are  discussed.  On  the  base  of  this  discussion  it  is 
made  an  effort  to  classify  the  gasdynamics  problems  in  which  clusters  play  role  and  to  propose  the  proper  quantum 
chemistry  method  for  each  class  of  problems. 


INTRODUCTION 


The  most  important  data  needed  for  gas  dynamics  that  could  be  gained  with  the  help  of  quantum  chemistry  is  the 
information  about  vibration  and  rotational  energy  levels  as  well  as  the  depth  of  potential  energy  well.  In  our 
investigations  we  have  used  a  number  of  different  methods  which  can  be  in  general  classified  into  three  brunches. 
These  are: 

i)  ab-initio 

ii)  semi-empirical 

iii)  empirical  potential  methods. 


AB-INITIO  METHOD 


Ab-initio  method  is  the  most  precise  but  most  computer  power  and  time  consuming.  Hence  usually  it  is  used  for 
simulating  relatively  small  systems  but  with  high  resolution  of  energy  levels  and  deep  elaboration  of  potential 
surfaces. 

In  the  reaction  of  H2SO4  formation  from  SQ  were  found  new  channels  of  formation  of  sulfuric  acid,  providing 
that  there  are  water  dimers  available.  There  was  discovered  a  new  transition  state  (1),  involving  two  water 
molecules  and  SO3  molecule,  having  much  smaller  energy  than  transition  state  with  only  one  water  molecule  (2): 


SO3  +  {HiO\  ^  H2  SO  A  ^HiO^HiSOa  +  HiO  (1) 

Were  checked  four  possible  channels  leading  to  the  formation  of  the  sulfuric  acid  (see  the  Figure  1): 

SO2  +  H2O  ^S02^H20  ^H2SOa  (2) 

SO2  +  {Hp\  ^  H2 SOa  •H2O  (3) 

S02*H20  +H2O  -^H2S0a*H20  (4) 

S03*(H20)2  ^H2S0a*H20  (5) 
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with  these  reactions  having  activation  barriers,  shown  in  the  Table.  1. 

This  set  of  reactions  was  first  treated  by  Morokuma  and  Muguruma  [1]  and  recently  repeated  by  us  [2].  Thus  the 
most  efficiently  the  formation  of  sulfuric  acid  would  proceed  with  the  help  of  reaction  (3),  which  assumes  that  we 
have  as  the  initial  reagents  water  dimers  and  SQ. 


TABLE  1  Activation  Barriers  for  reaction  (2-5) 

Reaction 

(2) 

(3) 

(4) 

(5) 

Energy,  Kcai/Mole 

32.2 

0.7 

5.3 

12.4 

All  the  calculations  were  made  with  the  help  of  the  program  GAMESS  [3], 


SEMI-EMPIRICAL  METHOD 


The  semi-empirical  method  uses  some  simplifications,  which  allow  treating  somewhat  larger  systems  and  fitting 
vibrational  properties  of  simple  systems  to  known  spectroscopic  data  (thereby  choosing  optimal  empirical 
parameters).  Semi-empirical  methods  are  especially  suitable  for  organic  systems  due  to  the  nature  of  C-H  bond.  In 
our  investigations  of  burning  problems  in  supersonic  aircraft  engine  and  flows  of  methane  in  industrial  conditions 
we  have  modeled  the  carbon  clusters  Q,  with  n  up  to  12.  It  appeared  to  be  two  main  types  of  isomers:  linear  and 
cyclic,  with  under  n=9  the  later  having  smaller  energy.  As  to  gasdynamics  consequences  it  was  found  that  in 
industrial  conditions  a  number  of  quasi-staedy  states  for  linear  and  cyclic  clusters  exist  before  the  steady  state  of 
cluster  distribution  on  size  is  reached.  It  is  important  to  take  into  account  for  this  effect  while  dealing  with 
technological  processes  where  carbon  clusters  take  part  (the  primary  stage  of  soot  formation  process,  the  initial  stage 
of  heavy  carbon  clusters  formation). 


Carbon  Linear  and  Cyclic  Ciusters  n=2  ...  12 


FIGURE  2.  The  energies  of  linear  and  cyclic  carbon  clusters  with  rFl2,  calculated  with  semi-empirical  method. 


EMPIRICAL  POTENTIAL  METHOD 


In  the  third  method  mentioned  -  empirical  potential  method  further  simplifications  are  made  through  the 
substituting  the  Schrodinger's  equation  by  Newton's  ones  and  introducing  to  these  equations  empirical  potential 
energy  already  containing  potential  energy  well.  Thereby  all  quantum  effects  are  supposed  to  be  plugged  into  this 
potential.  These  simplifications  allow  treating  rather  large  systems  (  up  to  10000  molecules  in  cluster)  and  what  is 
important  to  look  at  the  medium  region  of  cluster  sizes  corresponding  to  the  transition  between  clusters  and 
continuum  limits.  With  the  help  of  this  method  we  have  investigated  large  water  (up  to  100  molecules  in  clusters) 
and  large  methane  (up  to  100  molecules  in  clusters)  systems  [5-6].  It  should  be  mentioned  that  within  the  frame  of 
this  method  for  water  and  methane  clusters  the  dissociation  energy  obeys  the  linear  dependence  on  n  -  number  of 
molecules.  This  allows  us  to  extrapolate  its  behavior  to  at  least  next  few  n.  We  can  also  approximately  treat  the 
normal  modes  problem  within  the  frame  of  solid  state  physics  and  calculate  the  reaction  rate  constant  for  the 
reaction  of  adding  of  one  molecule  to  the  large  cluster.  Thereby  taking  into  account  the  opposite  reaction  of 
dissociation  with  the  help  of  gasdynamics  we  can  find  the  dynamical  distribution  of  small  droplets  and  aerosols  in 
the  flow  under  consideration. 


Methane  Clusters  n=3  ...  100 


FIGURE  3.  The  energies  of  methane  clusters  with  3=n=21,  calculated  with  empirical  potential  method. 


Water  Clusters  n=2  ...  100 
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FIGURE  4.  The  energies  of  water  clusters  with  2=n=21,  calculated  with  empirical  potential  method. 


All  the  calculations  ere  made  with  the  program  Tinker  [6]. 


CONCLUSIONS 


Summarizing  up  it  is  worthwhile  to  point  out  that  in  each  particular  problem  under  consideration  the  choice  of 
method  for  calculating  of  clusters  is  dictated  by  the  characteristic  sizes  of  clusters  as  well  as  how  detailed 
information  do  we  want  to  get.  For  large  clusters  the  linear  dependence  of  dissociation  energy  of  n  makes  it  valid  the 
application  of  empirical  potential  methods  and  molecular  dynamics.  For  small  clusters  both  ab-initio  and  semi 
empirical  methods  can  be  applied  with  the  former  giving  as  a  rule  more  precise  results  and  allowing  even  further 
corrections  with  the  help  of  Moller-Plesset  perturbation  theory.  For  organic  containing  systems  semi-empirical 
methods  are  the  most  convenient. 
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